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Program Summary 

Title of program: POLRAD 
Version: 2.0 April 1997 
Catalogue identifier: 

Program obtainable from: on request from e-mail; aku@liep.by 

Computer for which the program is designed and others on which it has been tested: 

Computers: all 

Operating systems or monitors under which the program has been tested: all 

Programming language used: FORTRAN 77 

Memory required to execute with typical data: 1MB 

No. of bits in a word: 32 

No. of processors used: 1 

No. of bytes in distributed program, including test data, etc.: 300 kB 
Distribution format: default ASCII else uuencoded compressed tar file 
Other programms called: 

PATCHY [1] — part of CERNLIB 
MINUIT [2] — part of CERNLIB 

Keywords: polarized particles, inclusive and semi-inclusive deep inelastic scattering, QED and 
electroweak radiative corrections, structure functions, higher order corrections, experimental 
data processing. 

Nature of physical problem: 

First and higher order QED and electroweak radiative corrections to the inclusive and semi- 
inclusive polarized deep inelastic scattering; experimental data processing. 

Method of solution: Numerical integration of analytical formulae. 

Restrictions on complexity of the problem: 

Only selective experimental cuts are possible. For (P(q;^) order correction only leading contri¬ 
bution is calculated. Electroweak correction is calculated for longitudinally polarized target. 
Typical running time: 

The running time depends on the options used. For example: 1) calculation of the total QED 
correction takes about 4 seconds of the CPU time per one kinematical point; 2) calculation of 
electroweak-|- (P(Q;^)-|-model for (72 7 ^ 0 takes up to 300 seconds per one kinematical point. 
References: 

[1] H.J.Klein, J.Zoll, PATCHY Reference Manual, March 1988. 

[2] F. James, MINUIT Reference Manual, March 1994. 
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1 Introduction 


Data processing of the modern experiments on deep inelastic scattering (DIS) of polarized 
leptons on polarized nuclear target requires correct account of the radiative corrections (RC). 
Our program POLRAD 2.0 based theoretically on the original approach proposed in ref.|^ and 
developed in the ref. was created to suit the demands of the present and future experiments 
with hxed polarized nuclear targets and at collider. Along with the possibilities of the previous 
versions of POLRAD |Q], which calculated the QED lowest order RC to DIS of polarized 
leptons by polarized nuclei, the current version gives an opportunity to take into account 
both electroweak and higher order effects and to calculate the RC for semi-inclusive polarized 
experiments. 

In section^ we present the detailed description of the theoretical basis of POLRAD 2.0 along 
with the explicit formulae. We start with the calculation of born cross section in subsection 
|2.1.1| . In subsections p.l.2| and |2.1.3| we present the review of the basic formulae for each of the 
radiative tails: elastic, quasielastic and inelastic, — and consider the case of ultrarelativistic 
approximation to the lowest order QED correction. The contribution of o? order correction is 
calculated in subsection |2.1.4| on the basis of structure functions formalism. The expressions for 
one-loop electroweak correction within the framework of standard theory and QCD-improved 
parton model are given in subsection |2.2|. The POLRAD 2.0 part that calculates RC in semi- 


inclusive case is the modihcation of the code SIRAD Iffl and is described in subsection 2.3 


The current data processing iteration procedure discussed in subsection p.4| is extended in 
comparison with the one in the previous version of POLRAD to £t the data using the CERNLIB 
package MINUIT. 

Appendix A is devoted to structure function dehnition, parameterization and models used. 
Most cumbersome parts of explicit formulae are presented in Appendices B and C. 

The current version of POLRAD gives the opportunity to choose one of hve nuclear tar¬ 
gets along with the type of its polarization, type of scattered longitudinally polarized charged 
leptons. It is possible to operate in standard SMC [Q, HERMES 0, E142 kinematics or 
to choose any other kinematics. Procedure of RC of experimental data can be organized with 
implementation of iteration procedure. The program is realized on Standard FORTRAN 77 
and does not require any changes when used under different computer platforms and operation 
systems. 


2 Theory 

We consider the process of DIS of longitudinally polarized charged leptons on longitudinally 
and transversely polarized nuclear target 

/ + AT r + X, (1) 

and semi-inclusive DIS (SIDIS) process when a hadron is measured in coincidence with the 
scattered lepton. The physical interpretation of the experimental data requires the separation 
of the Born cross section from background contributions known as radiative corrections, which 
originate from loop diagrams and from processes with the emission of additional real photons. 
Radiative events cannot completely be removed by experimental methods and so they have to 
be calculated theoretically and substracted from measured cross sections. 
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It is well-known that there are three scattering channels of virtual boson ( 7 , Z) on nucleus 
in dependence on transfer energy z/ = Ei — E 2 , where Ei{E 2 ) is initial (scattered) lepton 
energy: elastic, quasielastic and inelastic. Representative plot of dependence of the scattering 
cross section on u and square of transfer momentum = —q^ is shown on hg.l (only the 
regions that give sufficient contribution to RC calculation). Peak in the range I (for z/ = 0, if 
to neglect nuclear recoil) corresponds to elastic scattering. In this case the nucleus remains in 
the ground state. Range II stands for the quasielastic scattering i.e. direct collisions of leptons 
with nucleons inside nucleus. Wide maximum in the energy spectrum originates from the own 
movement of nucleons. Range III of inelastic scattering occurs when transfer energy is greater 
then pion threshold. 

On the born level both z/ and are hxed by the measurements of the scattered photon 
momentum. Hence, the channel of scattering is hxed too. However, on a level of RC radiated 
real photon momentum is indehnite, hence, v and are arbitrary so each of three channels 
contributes to cross section. Integration over the photon phase space may be presented to that 
in a plane of z/ and Q^. Adding the virtual photon contribution cr^, we have for the RC cross 
section 


a = a*" + a"' + + ah (2) 

Here each a denotes the double differential cross section (Ea/dxdy, and x, y are usual scaling 
nucleon variables, a®", a®^, are contributions of radiative tails from continuous spectrum 
(IRT), of the elastic scattering radiative tail (ERT), of the radiative tail from the quasielastic 
scattering (QRT) respectively. Also the contribution of electroweak correction calculated in the 
quark-parton model is contained in a®”. Both a and corrections are taken into account in 
a^n,ei,q,v^ To Separate the contributions we introduce the lower index, f.e. a®®® = a}®® -|- a^®®. 

All above mentioned contributions are valid in the case of inclusive scattering. However, in 
semi-inclusive case the transfer energy is above pion threshold, so in RC calculation one have 
to take into account only a®®® and a®®. 


2.1 QED formulae for inclusive case 

2.1.1 Born contribution 


Using (|A.1|) for the DIS cross section on the Born level, we obtain 


da SSr. 


dxdy Xs 


(Q'^ - 2m2)3i + (SV - 


^3 


'\ 2 \ ... ^4 


EmMEi (^2(Q ^7 - qy k20jp + Q J 

+(Q'^ - 2m^)(Q'^ - + (SX - 


— - ( -I- 4m^ -I- 12 rjki 7^2 ) ^7 “ x(-^ + S rik 2 )qri 


^8 

M 2 


( 3 ) 


Here ki{k 2 ),^,m are initial (hnal) lepton momentum, polarization vector and its mass respec¬ 
tively. Invariants are dehned in a standard way 
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( 4 ) 


S = 2kip, X = 2k2P = (1 - y)S, = -(fci - k^)^ = xyS, 


S^ = S-X, Ql = Q^ + 2m^, Sp = S + X, Xs = - 4m^M\ 

Pl is initial lepton polarization degree. An explicit form of hadronic tensor and generalized 
structure functions Aj are presented in Appendix |A.1| . 

The equality (^) is true for the any direction of polarization vectors and is exact: no ap¬ 
proximations were made yet. 

The 4-vector of target polarization p is the covariant representation of target polarization 
vector n. In the lab. frame p = (n,0), and n can be expanded in three components: parallel to 
initial lepton momentum ki - Ul, normal to ki in scattering plane (fci, /C 2 ) - and normal to 
scattering plane - n±. If a target is polarized along then we speak about the longitudinal 

(transverse) polarization. In the third case we speak about the polarization being normal to 
scattering plane. 

Let us build a basis in the 4-dimensional space. The process of inclusive DIS is determined 
by three vectors of incoming (outgoing) lepton ki{k 2 ) and of incoming nuclei p defining a 
hyperplane in the 4-dimensional space. We can choice the orthonormal basis system in this 
hyperplane {j^,Pl,Vt), where 

Vl = A;^/2(2Mfci - fp), 

^ {-SX + 2M^Ql))ki + Xsk2 - {SQ^ + 2m‘^S,)p 

Xl^\SXQ^ 

Among all possible basis system, our system differs from others in following: two space-like 
vectors pi and pt in lab. frame have the form (fi/^t,0), where are above-considered 3- 
dimensional polarization vectors. The basis vector system is uniquely fixed by this requirement. 
Basis in the 4-dimensional space is produced by adding to the system a 4-momentum p±{{n±, 0) 
in the lab system) orthonormal to the hyperplane. 

As a result for any 4-vector p we have expansion 


p 


pp p 

MM 


(mL)’lL - {Vk)’h - 


( 6 ) 


If p is target polarization vector, then for three above-mentioned cases we hnd: p = pi 
(longitudinal polarized target), p = pt (transversely polarized target), p = p± (target polarized 
orthogonal to the scattering plane). 

Initial lepton is always longitudinally polarized (for experiments considered). Using the 
expansion (^ for this vector we obtain 


^ = ^L = A//2(Afc^_2mp). (7) 

We note here, that calculation of real photon contribution requires to integrate over d^k {k - 
real photon momentum) some expressions, containing the scalar products of k and polarization 
vectors ^ and p. Since scalar products k^, kpi, kpt are easily expressed in terms of invariants, 
then our treatment allows to eliminate the intricate and tedious procedure of tensor integration 
used in ref. and signihcantly simplihes results. This is the most important advantage of the 
considered treatment. 
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By applying the ultrarelativistic approximation 


< S,X,Q^ 

and making transfer to scaling variables x and y we further hnd 

^ + (fi - ^ 62)(1 - 9) 

-FLFjvxy( 2 -y)yi) 

in the case of longitudinal polarized target and 

da dvra^S' 2 , ^ 

x\/xy(l — y)M 

- 2 P,P^ V (i/gi + 292)) 


( 8 ) 


( 9 ) 


( 10 ) 


for transverse one. 

2.1.2 Exact formulae of the lowest order 

The model independent RC of the lowest order can be written as the sum of bremsstrahlung 
and loop effects: 

a = a"^ + af + af + a'^. (11) 

The explicit form for these contributions was obtained in ref. . For the infrared free sum 
of a^ and we have 


a\ + — —d^ao + a'^p. — —{5^^ + d^ert + + ^vac)^o + • 

TT TT 


( 12 ) 


ap is the infrared free part of the IRT cross section 

Fmax Q ( Rmax j y-, 

aP = -a^y [ drY^loair) [ — 


2=1 


0 


55,(0, 0 ) 


R, 


/ dR 


(Q 2 + Rrp Q4 

Jp-2 


(13) 


(Q 2 + Rrp 


%{R,t) 


The integration region on variables R = 2pk and r = k{ki — k 2 )/pk is sketched on hg.|^. The 
limits of integration are dehned as 


R'm.nT. 


- (M + m^p 

1 + r 

92 I A T\j2r\2 


) ^max,min 


‘S'x ± y/Ag 

2 iW^ ’ 


( 14 ) 


Aq = ^^ + 4 MV, w‘^ = s^-q^ + mP 

where m,, is the pion mass. The explicit form of functions (r) is given in [Appendix B 
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The quantity appears when the infrared divergence is extracted in accordance with the 
Bardin and Shumeiko method [§ from The virtual photon contribution consists of the 
lepton vertex correction d^ert and the vacuum polarization by leptons and by hadrons 
0. These corrections are given by formulae (20-25) of ref.[^. Here we give ultrarelativestic 
formulae (m —0) for sum of and 


s:IR 


, X 3, ^ 1,2^ SX-Q^M^ 

+ Overt — Oinf 2 ^ In ^ T kl2 g'X' 


71 

T’ 


(15) 


where S' = X + X' = S — Q^, Im = InQ^/m^ and LB is Spence function (dilogarithm) and 


2\2 


^inf (^m 1) In 


S'X' 


(16) 


In the case of elastic scattering the nucleus remains in the ground state, so we have an 
additional relation 

R = R,i = {S^a-Q^)/{1 + ta) (17) 

resulting in 


e/ 


(Pa''’- 


ay 


'^At) 


8 ki 


2MlRir^ 


AdxAdy 


'^An 




^f{Rel,TA). (18) 


Here invariants with the index ”A” contain the nucleus momentum pa instead of p {p\ = M\, 
Ma is nucleus mass). The quantities 'Af are given in Appendix A-l| . 

Quasielastic scattering corresponds to direct collisions of leptons with nucleons inside nu¬ 
cleus. Due to self movement of nucleons we have no additional relation like (0)- As a result 
we have to integrate numerically both over R and r 


al = 


a^y 

' A 


Rl 


J / 


dR 


W 


i -2 


i=lj=l 


R"^ 

min 


(Q2 + Rr)2 


A?(i?,r). 


(19) 


The quantities Af can be obtained in the terms of quasielastic structure functions (so-called 
response functions, see Appendix |A.1| for explicit result), which have a form of the peak for 
uj = Q‘^/2M. Due to the absence of enough experimental data the fact is normally used for 
construction of the peak type approximation. The factors at response functions are estimated 
at the peak, and subsequent integration of response functions leads to results in terms of 
suppression factors Se,m,em (or of sum rules for electron-nucleus scattering [pT[| ): 


at = 


ay 


^ 4 ki 2M‘^ R^~‘^ 




( 20 ) 


To take into account the effect of radiation of many soft photons a special procedure of 
exponentiation was applied |^. In the code it is realized by the following substitutions: 


a 


el 


'y‘^{l — x/Af 


1 — xyj A 


a 


el 

1 5 




'y\i 


X 


xy 


af, 


(1 

TT 


exp 


^infj (1 

/ TT 


( 21 ) 


a, 

where Q = — (/^ 

TT 


!)• 
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2.1.3 Ultrarelativistic approximation 

To simplify and accelerate the procedure of experimental data processing when rapid analysis 
is more important than accuracy, it is convenient to have the approximate formulae. In the 
case of RC calculation one can choose ultrarelativistic approximation: 

m^,M^ <^S,X,Q^, ( 22 ) 

that allows to calculate exactly hrst two terms (corrections ~ aim and ~ a) of expansion over 
the leptonic mass m of the lowest order cross section. 

The expressions derived under such an approximation are compact and have good accuracy. 
Also this approach allows to avoid numerical uncertainties when the value results from the 
difference between two large and sometimes inhnite quantities, that is especially signihcant 
when quadruple polarization is considered. 


A. Inelastic radiative tail 


Considering r-dependence of quantities %(r) in ( p!3D one can see its peaking structure, so called 
s- and p-peaks H (or ki— and fc 2 ~peaks if follow [^): O^j (r) ~ 9*,(r) + (r). Using the 

identities 

9yW3i(-R.T) = 9*,(T)(9i(fl.T) - + 9*,(T)9i(fl,T,). 

(23) 

ef,{TmR,T) = «,yT)(3i(fi,r) - + «pr)3i(fi,rp), 

one can extract and analytically integrate the terms corresponding the mass singularities. The 
hrst terms in right-hand sides of (|23| ) are free from mass singularities and so they contribute 
only to ~ a correction. So one can adopt = 0 before the integration over r (or 

photon radiation angles). SF’s do not depend on r in the rest terms. Hence, the last ones 
can be integrated analytically using the methods [^| and contribute to the leading correction 
~ aim- 

Results for the infrared free sum of contributions from inelastic radiative tail and loop effects 
could be presented as following 


a 


rn I V _ " Xin ^ , V , ^in , in , „ 

<JI + CTi — —0i ctq T CTj, <Jg ajj a. 


TT 


Factorizing terms of the total correction have the form 


^i" = 2 [(4 


- ln“(l -y)- 1|, 


where 


(24) 

(25) 

(26) 


Ssp ^ 21ii((l - 2p)(l - Zs)) + 3, 

The quantity is correction due to vacuum polarization effects by leptons and hadrons: 

T — iKaz + ^vac)^0- (27) 

Quantities cr**), are the contributions including the second terms in right-hand side of 

1 


^s,p = 


a 


dz 

1-z 


1 + z'^)lm - ‘2z]as,p - 2{lm - l)cro) , 


(28) 
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The quantities as,p are expressed in terms of born cross section ao = cro(S', X, Q"^): 

as = yao{zS,X,zQ‘^)/{z-l + y), ap = yao{S, X/z^Q"^ / z)/z{z - I + y). (29) 


The low limits of integration in (|^) are 

- ^~y 

1 — xy 


, Zp=l-y + xy. 


(30) 


The mass singularity free terms in (^) contribute to cr*". After splitting the unpolarized, 
polarized and quadrupolarized parts, we have 


^in 

= 


2ay / -I {Tt + PLPN{T^^^ + T;^^) + 

X ^ 

in the case of longitudinally polarized target and 

< = 2ay I p {Ti + PlPn(.T^\ + T/j.) 

X ^ 

for transversely polarized target. Quantities T are 


(31) 


Qn '~pL 
6 


(32) 


Tt = A(Fi) + ^A(F2), 

^pII ~ 

rpx _ -\- Ux)gi{^itx) 

^11 “ ’ 


T,^ = A(ti) + ^£ 2 ( 62 ), 

+ 2 A(S 2 )/Q|), 

4^Mm(Q|^i(^, 4) + 2^^2(C, tx)) 

~ XQ|(^XQ2)1/2 


where 

C,{X) = 2L^ - + Li + - Li + 2L„ 

C,{X) = {T- QlSx)L^ - {QlX + T)Li + (T - Q|^)Li, 
/:3(^) = 2SpL^ - {us + X)Li - QlL[ + {S + u.)L’; - QIL^ 
A(^) = -4:XL^ + (2X - 5)Li - n,Li + (2Q2 _ ^Ql^ 


— {X + 2ux)L^ — UxLI'x ~ '^xLx + 2SxLt, 


A(^) = {SxQl - Sl)L^ + {SSp + 2XnQLi + n^XL^T 


(33) 


(34) 


-( 2 ^n, + SpX)Ll - UxSLl. 

Here u = S^ — <5|, = S — Q|, Ug = X + Q|, Q| = T = S'^ + X^. Quantities L, dehned 

below, have to be computed for the same argument T {T = F 12 , ( 71 , 2 , & 1 - 4 ) as in C-{T\ 
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m,t) 


jl ^limits,x) ^s,x 

^S,X ~ n.D2 


t2 


+ 


Vjq 


uQ^ u 


dt 


j^k _ d~i{^,ts,x) '^s,x 


t2 


s,x 


L^ = I dt 


+ 


S, X uQl S, X 


dt 


- ^ii^,'ts,x) 

\t-ts,x\ 

ts,xXi{i,t) - tXi{i,ts,x) 

tts^x\t tgx I 


( 35 ) 


t2 


11 = 


Ux{t - t^){Xi{i,t) - Xi{i,ts)) + Us{t - ts){Xi{i,t) - Xi{i,t^)) 


ts\ "1“ ^s||^ ^x 


dt 


m,t) 


Ua 


t2 Ug 


t2 

C " 
ts,x 

-v 2 


dt 


m,t) 

t2 ■ 


The integration region on variables ^ = —q^/2pq and t = —(f {q = ki — k — ^ 2 ) are plotted on 
fig.^. The limits of integration are 


ts,X Q 


2{M1_ 


u{Sx: ± sJXq) + 2M2g2 
2{u/^ + M2) 


t2,l — 


(36) 


B. Elastic and Quasielastic radiative tails 

In the case of ultrarelativistic approximation in calculation of RC from elastic radiative tail, 
due to strong dependence of formfactors Fj on the square of transfer momenta Q^, the leading 
contribution to the total cross section gives only t—peak (Compton peak) and contributions of 
s— and p— peaks are suppressed [Q. For af we have 


af = a^ + PlPnK + 


(37) 


where index A corresponds to the considered nuclei and m, p, q dehne the unpolarized, polarized 
and quadrupolarized contributions. For spin 1/2 nuclei Qn = 0 and for spin 0 nuclei Pn = 
Qn = 0. 

For proton, deuteron and carbon we obtain the results 


= 


aP = 


(T,. = 


= 


a 


IT 


dpA 


S J r]A 

'Hmin 


dpA 


cA 

—Y_ 

S J r]A 

Vmin 


(A'(Ff + ,,^f|) - (F, + F-S-), 


(Fi + F2)(xXF, - (Fi + F 2 )), 


a" 

'S 


00 

n f -^{(F! + lF^vl + lFlr,AX-Hl + VA)Fl], 


(38) 




dpA 


—Y 

S - J VA 

'Hmin 


Fm(|(l + r]A)Fm - {Pc + iFqPA + \FmPA)xX 
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<i 

- 


a" 

'S 




^ J 

nin 

XX 

- ^ Fq{3Fc + SrjAFm + rjAFg) 

(1 + Va) 

-27]AXFq{AFc - ^xFra + fr^^Fg)|, 


= 


a" 

~S 


Z^Y_ 


Vmi' 


drjA 

Va 


XF^ 


for the case of longitudinally polarized target, where 


Xi = X + 4xr]A - 4:r]A, X = 


2r]Ax‘^ 


Y+ = 


l±{l-yf 


y 


t 


VA = 


4M1’ 


^min — 


X 


(39) 


and formfactors 


Fi = 


Ge + VaGm 


Fq = 


4(1 — x) ’ 
Gm — Ge 


(40) 


I+ VA I+ VA 

Polarized contribution in the case of transversely polarized target is proportional to nuclear 
mass and therefore equals to zero in the case of ultrarelativistic approximation (^). Neverthe¬ 
less we present an explicit formulae for the hrst nonzero order correction: 
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a xy 


p _ « M 7 dr]A f 

Vfnin 


J Va I V ^ yVA^ 


X \ (2 — y)x 


y 


X + 


(41) 


where yi = 1/y — 1. The cross section dependence on nuclei formfactors are contained in 
quantities Q. For protons and deuterons we take 


gi = {F^ + F^)\ = F2(Fi + F2), 


Ql — ~x(l + Va)F^, ^2 — Fm{Fc + \VAFq — \Fm)- 


(42) 


Also we present relations between quadrupolirazation parts for the cases of longitudinally and 
transversely polarized targets: 

= -Kii = (43) 

For the calculation of the contribution to RC from the quasielastic radiative tail in the case 
of ultrarelativistic approximation we use the results for the proton target for elastic radiative 
tail replacing 




SemF^m, Fl{Q^) SM^)Fl{Q% 


Fe(Q')F^(Q' 


Sem{Q‘^)Fe{Q‘^)F^{Q^). 


(44) 
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2.1.4 Higher order effects 

There are no known reasons to consider the 0{a^) corrections to be negligible. We follow the 
method of structure functions used in for the calculation of the higher order electromagnetic 
radiative corrections to neutral current unpolarized lepton-proton DIS and generalize it for the 
case of polarized leptons and polarized nuclei target in the current version of POLRAD. For the 
case of s- and p-peaks the formulae could be obtained in the terms of the Born cross section and 
practically coincide with the expressions for unpolarized particles, however the contribution of 
the f-peak which is extremely important in the cases of elastic and quasielastic radiative tails 
has to be obtained for polarized DIS. 

The sum of second order inelastic correction and correction due to loop effects has the form 


An I £in ^0 , An , An , An , i i i i i 

<^2 +^2 = ^2 +(^Vs + (^Vp + ^ss+(^pp + (^sp + ^ls +(^lp +^fs + ^fp^ 

where factorized part is 

^ A(Q") + •iiraS.AQ") + I'lK + 4Li2(l - Z,) + 12Li2(l - Z,) - 


(45) 


(46) 


and5(Q2) + 

The contribution from vacuum polarization if coincide with real photon radiation has the 
form 


1 


dz 


^in _ ; I 

^vsyp-^A^J 1-z 

• 2 s,p 


((1 + z^)d(tx,s)o's,p — 2S{Q‘^)aoj], 


(47) 


where = zQ^ and tg = jz. 

Next three terms correspond to the cases when two radiated photons are collinear to incident 
electron (uss), outgoing electron (cTpp) 


in _ ^ j 2 

’^pp,ss — 


dz 


(21n(l - z){l - ZpA)) - ln ;2 + 3)((1 -1- zAp,s - 2ao) 


1-z 

+ ((1 + z)\nz — 2(1 — z))ap^ 


(48) 


or when one photon is radiated in incident electron direction and the other in the outgoing 
electron direction [a 


sp J 


in ^ Ap } 
^rn 


dzo 


1- Zi J 1- Z2 

2p(2l) 


(1 -|- Z2){1 + zl)(Tsp — 2(1 -|- zl)as — 2(1 -|- zl)ap -|- 4(To 


Here 


Zs(z) = (1 - y)/(z - xy), Zp(z) = (1 - y + zxy)/z, 


(49) 

(50) 


and as,p depend on 2 : 1^2 and are given by (^) with 2 : —> zi ^2 respectively. The quantity agp 
depends on both zi and Z 2 and also is calculated in terms of born cross section (Tq: 


(Zsp = yuoAiS, X/Z2, zA/Z2)/{z2{z2Zi - 1 + y)). 


(51) 
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There are two channels (singlet and non-singlet) of fermion pair prodnction that give a 
contribntion to order RC. The singlet channel 

2 1 

+ z)\nz + l-z + l{l- z^)z)ap^., (52) 

^p,s 

corresponds to the case when incident and ontgoing lepton as well as leptons of the nnregistered 
pair belong to different leptonic lines connected by additional virtnal photon. 

The rest non-singlet part 


a 


in 

fpjs 


127r2 




91 

m} 


I * 

Zp,S 


9 + z^) 

{l-z) 


^P,S 


(53) 


arises from the two-lepton decay of additional virtnal photon. 

The main contribntion to second order elastic and qnasielastic radiative tail arises when 
additionally radiated photon is collinear to incident or ontgoing fermion line: 


^el _ 7 X I ^eL , ^ei , ^ei 

^2 — -^^mOspC^l + CTy + 


a 


^el 


^el 


^el 


el 


^el 


^s,pt 2vr"* 


^^ {1 + z^)atlp - 2at 
^ l-z 


el 


(54) 


where the qnantities are obtained in terms of approximate elastic radiative tail (^71) af = 
af{x,y,S)-. 

= y(Tf{xyz/{z + y -l),{z + y - 1)/^, zS)/{z -1 + y), 


(55) 

= y(^ii.xy/{z + y-l),{z + y- l)/z, S)/z{z -l + y). 

The integral over ^ in (Q) can be calculated explicitly. For longitudinally polarized target the 
results have the form of eq.(|S^) 



(56) 


where the quantities are given in [Appendix Q . The functions quadratic com¬ 

binations of nuclear formfactors and could be found by comparison with (|HB|) which can be 
written in the common form 


a 


(J. 


u^p 


= —Y- 


S 


± 


Pmin 


dr]A 

VA 


A2 




(57) 


The correction due to vacuum polarization ay 
tional factor ^h(4M^?7^) under integral. 


is dehned by formulae (^) and (|38D with addi- 
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2.2 Electroweak radiative correction 


The next evident step both from the theoretical and experimental points of view is the treatment 
of electroweak effects contribution. So we included in POLRAD 2.0 the results of ref. IITO 


for one-loop electroweak correction within the framework of standard theory and the on-shell 
renormalization scheme in t’Hooft-Feynman gauge. The result for the correction is obtained as 
the sum of loop and radiative effects 


+ ^Vl + (TVq + (Tbox + ^ I] <^ 9^0 + 




(58) 


where af is the correction to boson propagator, avi,q, (Tbox are infrared free parts of lepton 
and quark vertex functions and box graphs. The loop correction is calculated on the basis of 
ref.[|^. The quantity cx^ is an infrared free part of the real photon emission cross section. 
The correction 6q is obtained after infrared divergence cancelation. It is factorized front of the 
born cross section on a quark cxo and is an analog of quantity in (^). For radiative effect 
the methods developed in ref. ^ are applied. The quantity EgCr 
leptonic , a}^), hadronic radiation and their interference 


can be derived in terms of 


= 




ij='y,Z 




)} 


The quantities al\ have the form of one-dimensional integrals over ^ 


(59) 


al^ = 


4 I e 




+ R 


'O 


nFfK)-r?i'Ff(x) 


(60) 


where are the combinations of kinematical invariants, RyA ^^*4 ^va functions of elec¬ 
troweak coupling constants and parton distributions ||^. The hat-quantities in (^^ are not 
large corrections arising from the non-leading terms of the expansion of polarization vectors 
(see 13) • 


2.2.1 Correction to leptonic current in QCD-improved model 

An implementation of QCD-improved parton model for the most important case of leptonic 


current correction requires an additional generalization of (^ |60D valid for the simple parton 
model and cannot be generated directly because an analytical integration over Q\ been already 
done. 

As a result correction takes the same form as eqn. 


a 


^in I _ ^in I ^ \ \ ^in \ ^in i ^in 

lept + (^llept — — —^1 Cr0 + (^r + <^5 + + <Jj. . 


TT 


(61) 


The quantities a*” and 5^"' are dehned by (p8| ) and (pS]). The Born contribution can be written 


as 


na 


(^0= E ^0 = E 204 

iA='y,Z 


X{Y+RtF!^{x, Q2) + Y.iriFi^x, Q^)}. 


(62) 


*J=7 

Here the following notation is introduced in (|A.1CI|) 
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The contribution of electroweak loops can be found in the form 


cr^ = SR'l), 


( 63 ) 


where 




= -2YVRl^-2 YV^xR± 


R2{—Q‘^ ■, Mz) + n/ (1 — Pl)—A3(— Mw)], 


dvr 


6Rf = 6 rI'^ = -2{ir + n ^){Rf + rI^) - n^^(p]7 + xRi^) 




(64) 




+ (1 - Pl){—yA2(-Q^, Mw) + + af -v'i)JXs{-Q‘^ ■, Mw)}]-, 

4st, 2st s,„ 


SRl^ = -2U^Rl^ - 2W^Rf 


a 


+^[{yrR± + yrR^^)A,{-Q\Mz) 


+(1 - PL){vf + af){—A 2 (-g^ Mw) - 3^A3(-Q^ Mw)}]. (65) 


Here Mzw are masses of Z and W bosons and 


Q2 ’ g2 + ^|’ Q2 • 


( 66 ) 


Quantities are dehned by formulae (A.2,3.17,B.2-5) of and A 2,3 by (B.4,B.6) of [|n 

In the electroweak case the term a!” has form 


di r 


<= T. -X 


i,j='r,z X 


e 


R^T)^ + RXiTi + PL{,X^"/f^ + XyM!^) 


(67) 


where 


T|^ = 


4?Qt 


|<3|£,(n-’) + 2£2(F«)], 


2^0 _ 


4eg| 


A(i"^') 


( 68 ) 


m(^2 ± ^2 


_ _ V _ xJ pij 

^ 2^UxXQj 


Functions Pi, 2,3 are dehned in (^). 
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2.2.2 Correction to hadronic current 


Exact formulae for correction to hadronic current could be obtained from when b = h and 
could be easily generated for the case of QCD-improved model (see also ref-HI^). 

Result for leading log approximation could be presented in standard form PD[]: 


H , (69) 


where 




2^ 0^ 
1 


-ftix) In- 




(1 — xYm^ 

+/T{^^h-A-)-T4^hw}}- 


(70) 


Sometimes hts for partonic distributions are constructed from the data extracted without taking 
into account the of hadronic current correction. Therefore if such a £t is used in the calculation 
the correction does not have to be taken into account. 


2.3 Semi-inclusive physics 

We calculate the radiative corrections to data of semi-inclusive polarized experiments when 
a hadron is detected in coincidence with the outgoing lepton. In this case the cross section 
depends additionally on variable 2 ; dehned as 


z = 


P1P2 
PiQ ’ 


(71) 


where pi, p 2 {p^ = rn\) and q are 4-momenta of initial nucleus, coincident hadron and virtual 
photon. This variable corresponds to the amount of virtual photon energy transmitted to 
measured hadron in lab frame. 

For the Born cross section of semi-inclusive DIS we use the formula 




27ra‘^ 


CTn = 


dxdydz Sxy 


KE+(I, z) + PlPnFS1:-(x, z)] 


(72) 


where 

F^ = 2{l/y-l-yMx)+y, FS = y-2, = M^S (73) 

and z) are defined in |A.2| . 

The lowest order QED correction was calculated in ref. |^, ^ and can be written as the sum 
of factorizing and non-factorizing parts 


_ d^CJEM F I ® c s , V 


(74) 


where 

= by + - - In^ ^ -Im -2 + Li2(r2/ri)- 

^inf = (Im - tu = min{tim,hm/r+} 


(75) 
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and 


tim = 2/(1 - x) + hn - (M + TUhf/S, r± = [y + 2/1^ ± {y^ + 4/iAra;2/)^/^]/2r9, 
t2m = y{^- z), ri = 1 - x^y^ - yr^, r2 = - yrg, rg = 1 -//(I - x), 

ri = l-xy, rg = 1 - //, rg = /xa^x, vg = + y - xy. 


The finite part of ([72D has the form 

Q ^lu t'2u O 

v'5 n r 

I dti I dt2\ 

zxtyt 


= 


2a^ 

~S' 


^lu ^2u 2 

r dti [ dt2{^[FiT.-^{irz) + PLPNFf,T.-{xrz)]- 

J J Xf Vi 


hd i-2d 

Oitii - h) 
xy 


[F7^S+(x,z) + PAPivF/^S-(x,z)]}, 


where the integration region is defined as (see fig.|^) 

t\d 0 ) tlu ^XliTl{tlmit2m/^—'} 1 

t2d = tir_, t2u = min{t2m, tir+) 
and we select the nnpolarized and polarized F^ terms: 

Fr = {xt/yt){[-2yl‘^ - + X2/J](r5(l - ta) - ^NXt + yl/F)F 

+ [-2yP + + xyl]{r5 + fa - ^NXt + yljF)- 


-(1 + yr^)l^ + {rl + //rg)/' - 2/r^/°} + -xy{I^ - /') + l\ 


and 


= {\ly^{G\xyyt - 2xt] - yyt{t\ - t2)P + xV^2/+ 

+ [xy{y + 2 /t - 1) + (fi - t 2 ){yt /‘2 - 1)]/^- 

-[xy{yt - 1) + (fi - fa)( 2/^/2 - 1)]/'}. 

In displayed above formnlae 

y = vri^lS, j/+ = // + f2, Xi = X 2 / + fi-fa, yt = y - t2, 


^ Q ^ _^(j2 ^ j2^ ^ ^yj 


Here we nse the formnlae 

/» = /» = A’''^ /’ = [A'‘ + P = [B^ + 

P = !/[(l + 2r8)ii - (r, + xre}t2](P f, 

P = !/[(f5 - 2r8)(i - (rsr, - 


(76) 


(77) 


(78) 


(79) 


(80) 


(81) 


( 82 ) 
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where 

and 


A = ti- r4t2, B = r^ti - r2.t2, C = 4/i{r6tit2 - - rgtl}, 


rQ = y + 2fiN, 


rj = 1 — 2x. 


(83) 


(84) 


To take into acconnt the higher order contribntion of soft photons the standard exponentia¬ 
tion procednre is nsed 


SeM —>■ Sexp = eXp[(a/7r)5inf](5EM, 


so that 


a 


" = crj(l + i5„p). 


(85) 


( 86 ) 


Since the experimental analysis is often performed for all pions with z > zg we have to 
integrate the expression (^) over 2 ;. As 2 ;—dependence is contained only in qnantities and 
integration limits ( 0 ). then one can use the following identity 

[ dz [ dtidt 2 '^(pqDf{z) = [ dtidt 2 '^ipqDf{zo,yt/y), 


where 


Dq{zoX)=Cl dzD^{z), 


H/ 


'20/C 


(87) 


( 88 ) 


and ipq = ipq{S,x,y,ti,t 2 ) is an arbitrary function. Hence, substituting some new effective 
fragmentation function {zo,yt/y) instead of D^[z) and D^{zq, 1) instead of D^{z) one can 
represent the integration over in in the form of eq. (^) and get 


a^{x,y) = f dza^ = cro(l d^xp) \ ^ 

Bo 


Dq{z) —^ D^{zo,yt/y) 


D^iz) 


D^ZgX) 


(89) 


2.3.1 Experimental cuts 

In the previous section the integrated in the whole kinematical region over hadron variables (pn 
and p± cross section was considered. For the real situation the region of integration is limited 
by the experimental cuts. To deal with the ones on the angles of registered hadron a special 
procedure was developed. For cross section we have 


da 


R 


dxdydz 


1 d^k dz 
2 tt 2kn dz 


dp±d(()HO(,silX id max “ siu^ ld)9{sm^ -d - sin^ ^min) 


da 


dxdydzdp^dcpH 


, (90) 


where'd is the angle between the beam direction (ki) and hadron momentum in the lab frame. 

are removed. We note, that (pOD is valid 


One can obtain ([77|) straightway if d—functions in 
when asimutal symmetry of detector is supposed. As the distribution on p± is unknown it was 
approximated by d-function 

D{z,p±) = D{z)5{p±)^, (91) 

according to the normalization condition 

[ D{z,p±)dp±d(j)H = D{z)- (92) 
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The presence of <5—function in (PID takes off the integration over p±, so the integration over (pn 
becomes trivial as arguments of 6*—functions are now not dependent on p 2 (and on (j)H), but 
acquires dependence on photon momenta. 

Implying peaking approximation in (pOD we preserve only ti dependence for 6*—functions, 
that allow to carry out analytical integration over one of photon variables. As a result we have 
the same formula as (PD|), but we replace 


jl,2 jl,2 _ 


rsiy - d) 


l/« |9(i. _ r,{y - d)) 


(93) 


where we take into account, that in nonradiative case the approximation (|^) is equivalent to 
the assumption, that only those hadrons are registered, for which 


As 

fdmin A V A (djnaxi 

yl 


where 


A 


sin^ d. 


min,max 


mm,max 


2/ijv 


d = 


Aq 


1/2 




(94) 


(95) 


Ag = |/(-4a;/ijv/i - xy - py + x - pN^'^y)- 


2.4 Iteration procedure of data processing 

From the POLRAD beginning particular emphasis has been placed on the procedure of RC of 
experimental data. In the current version the iteration procedure, which allows to extract Born 
data sets for cross sections, SF or asymmetries from observables ones taking into account the 
radiative effects is realized both for the cases of inclusive experiments. 

As an example we consider the procedure of radiative correction to extract the structure 
function gi{x) from measured spin asymmetry A™ with error e*. The spin average structure 
functions are considered to be constant and g 2 {x) equals to 0. The measured asymmetry is 
dehned as 

yir = ||- + AA(9.), (96) 

where the radiative correction to asymmetry AAi can be written in terms of spin-average and 
spin-dependent parts (a"’^) of cross sections (^ 


AAi = 


a 




yin 

'u 


+<+<') 


a 


'((1 + KK, + < + rf + ai) 


(97) 


where 5^ = cr^/aQ = ct^/ctq. The Born and inelastic radiative tail polarized parts of cross 
sections depend on SF gi, and in the last case the dependence is non-trivial. So the equation 
(p6|) becomes functional one in gi. This functional equation transforms into a system considering 
the extraction of gi in concrete binning over x in n kinematical points Xi (i = 1, 


A 


m 

li 


dli 

Fi 


+ AAi(5fij; j 


l,...,n). 


(98) 
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Usually the iteration methods are used to solve such a system of equations. The variant of 
iteration formula is ambiguous, but in practice only two types are used. The first and most 
evident one is to take for k-th. step: 

= FMu - j = 1, ...,n)). (99) 


Another possibility to obtain the formulae for iteration procednre arises when both born and 
radiative correction cross section are separated into spin-averaged and spin-dependent parts. 
Then for the measnred asymmetry we have 


A™ = 


1 (T^ + af 


Dak 


0 + 


gi/Fi + al/Dal 
1 + ak/ak 


Thus we obtain the iteration formulae 


(k) 

9ii 


= F, 


Am 


l + % 


Un 




Da}, 


( 100 ) 


( 101 ) 


where in right-hand side the dependence on gi is contained only on the level of RC, bnt not on 
the Born level. 

ik) 

On the each step of iteration procednre gh is htted with the help of CERNLIB package 
MINUIT taking into account an experimental nncertainty €i{gi) = ei{A^)Fi. On the first step 
we adopt gf^ = The procednre converges within 4-5 steps. As a resnlt we extract valnes 

of gii and parameters of its fit. 

Besides, package MINUIT is used to £t the data with the account of experimental nncer- 
tainties that gives the opportnnity to theoretically calculate the error propagation of statistical 
nncertainty of fitted experimental data to the valne of AA: 




1 1/2 


6^(1 -F 5iu — 5ipf + ^ 51 






r _ ^in + + ^el ^ ^a^ dpk 

< + < hsp.dAf 


( 102 ) 




^iti. 




9(Tfn dpk 


9pk <9A; 


The snm runs over parameters pk of fitting fnnction f{x). The first derivative is calculated 
by direct calculation of a^^ by POLRAD with using instead of model Ai (or gi). The 

second derivative is obtained from the system of linear eqnations 


l\p 
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dpk 

(9A™ 


N ^ 

u? 

i=l 




" I 

fcn' 


X5 





(103) 


Let us present an explicit formulae for three most important cases, when proton structnre 
fnnction gi is extracted from the data with hydrogen target and neutron strncture function 
gi is extracted from the data with denteron and ^He targets. Also in the last subsection we 
consider the other target possibilities withont taking into acconnt experimental nncertainties 
(not using CERNLIB package MINUIT). 
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2.4.1 Proton, Deuteron and Helium-3 targets 

The formulae of previous section could be applied directly for experimental data on hydrogen 
and deuteron. Three-parameter £t is used for htting of the spin asymmetry |P: 


AP{x) = A + x^{l-e-^^), 
Afix) = {e-^^ - - x^). 


(104) 


For proton asymmetry the one- and two parameter hts of ref. 122] are also used: 


Ai{x) = Ai{x) = Ax^. 


(105) 


For the neutron asymmetry extracted from the deuteron and helium-3 data the nuclear 
corrections have to be taken into account: 


AUx) = 


Af- (1 - fd{x))PpAP^{x) 


fd{x)Pn 

Dilution factors for both cases are given by the formulae 

1 „ . . 1 


fd{x) = 


Fi/F^ + l' 


fd{x) = 


2Fi/F^ + V 


(106) 


(107) 


respectively. The numbers Pp and Pn are effective nucleon polarizations in nucleus. For deuteron 
target Pp = Pn = 0.5 — 0.75uJd, where ujd is the D-state probability (~5%), and Pp=-0.028 and 
F„=0.86 |2^ for hehum-3. 

In both cases we construct the £t for neutron asymmetry using Schaefer’s model (see 
Appendix |A.4|) . 


2.4.2 Other schemes 

In many experimental cases when the model for structure function extracted is unknown it is 
convenient to use simple spline approximation of experimental data. Below we consider two 
examples, hrst when gi and g 2 are extracted simultaneously from data on spin asymmetries Ai, 
A 2 and the second, when bi is extracted from the data on quadruple asymmetry Aq. 

A™ = Af (x) -I- AAi{gi,g2) = -—^ + AAi{gi, g2), 

. ‘ , (108) 
— /if (a^) + AA2{gi, g2) — —pzp- h A/l2(ili.<?2). 

_ 2 

where 7 = ^/Q jv. SF gi{x) and g 2 {x) are calculated on the each step: 

sh = TTT^ K - A/1.(9<"-0 ,j 70) +7(/l? - A/12 (j!"- 0 (»-i) 

(1 + 7^) V 

2' K - A/l2(ji"-0,!,7‘') - J(AT - AA^iat'Kst 

7(1 + 7 ^) ^ 
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For quadruple deuteron spin-dependent SF hi{x) we have: 



(110) 

and 



(111) 


Also the spline method, described above can be used for the cases of previous section to 
obtain the model independent data. 


3 User manual 

3.1 Program structure 

In this section we present the common structure of POLRAD 2.0 along with the short descrip¬ 
tion of subroutines used. 


3.1.1 Main program body 

The code POLRAD 2.0 operates under PATCHY from CERNLIB |^. It means that there are 
two hies POLRAD20.CAR and POLRAD20.CRA. The second hie contains a set of switches 
for user. The hrst one is text of code for all possible combinations of the switches. FORTRAN 
code is obtained by calling like pat polrad20, where pat is a simple executable hie. For examle 
for SUN station it can have form: 


ytobin $1 $1 $1 $1 .go 
ypatchy $1 $1 $1 $1 .go 
rm y.lis 
rm $l.pam 
rm $l.lis 

POLRAD 2.0 allows user to conduct the calculations and to obtain the results for the any set 
of the the following positions (in brackets we refer to the corresponding theoretical description). 

1) Exact calculation of the total RC to DIS of polarized particles (see subsection E3); 

2) Electroweak RC calculation along with 1) (see subsection p.2|) ; 

3) The contribution of order correction calculation along with 1) (see subsection 2.1.4 ); 

4) Approximate calculation of the total RC (see subsection O; 

5) Exact calculation of the total RC to SIDIS of polarized particles. 

Hence, the main hie POLRAD20.CAR consists from the parts (patches) POLRAD; EX¬ 
ACT; POLRAD_ADD; SIRAD corresponding an above mentioned positions, parts (patches) 
STRFFUN; FITS2; INTEGRAT of the common use. The correlation between the patches is 
presented on the hgs.H-^ 

Patch POLRAD. POLRAD guides the calculation and sends the obtained data in output 
hies. On the hrst step of program run subroutine^] CONKIN sets up the value of the invariants 
in dependence on the given lepton, target and kinematics. Subroutines FSPENS and FSPEN 

^Usually one subroutine corresponds to one DECK from PATCHY with the same name. 
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calculate Spence function. Subroutine DELTAS calculates the factorizing part of virtual and 
real leptonic bremsstrahlung. Subroutine BORNIN calculates the Born cross section on the 
basis of known formulae . 

Patch EXACT. In TAILS and FFU the r-dependent part of an integrand (%(t) in (1)) is 
calculated. The integration procedure has the following steps: 

1) QQT and QQTINT set up the integration parameters and call the integrators; 

2) Integrators use the integrands which are calculated in subroutines RV2DI, RV2 and 
PODINL. 


Patch SIRAD deals with semi-inclusive RC. Main program MAINPRG guides the calculation 
of Asymmetry or r{z). Subroutine COM VAR dehnes some often used variables and the limits 
for the phase space of the emitted photon. FSP and FSPl give the Spence function. Subroutines 
EXHH, QXT, INTEG and Nil carry out integration over photon variables. Models for parton 
distributions and fragmentation functions are given in DZ, FITDZ, SIGMA, QS QP. Subroutine 
DOS and DOP dehnes unpolarized and polarized Born cross section. 

Functions FYS(FOYS) and FYP(FOYP) dehne the sum and difference of cross sections with 
the target polarization vector parallel and antiparallel to the momentum vector of the incident 
lepton. Radiative effects are (are not) included. 

Functions FCS(FCP) and FRS(FRP) are the integrands of the unpolarized (polarized) part 
of bremsstrahlung cross section over invariant variable corresponding to the energy and to the 
polar angle of the emitted photon. 

Subroutines VPQRK and Z0_EXH calculates effects contribution of the quarks to the vac¬ 
uum polarization and effects of Zo-exchange at the Born level. 

Subroutines INPUT and OUTPUT create input arrays of kinematics variables and form 
output hies. 

Patch STRFFUN. Subroutine STRF carries out the calculation of inelastic, elastic quasielas¬ 
tic SF (Ai(R, r)), dehned for diherent cases in [Appendix A| as the combinations of usual unpo¬ 
larized (F 2 , Fi or R), polarized {gi, ( 72 ), quadruple (&1-4) structure functions and formfactors. 

SF F 2 is calculated in dependence on the user dehned model for the whole kinematical range 
in the subroutines: 

GOMFST and F2BRAS calculate F 2 outside and inside resonance region on the basis of 
experimental data from refs. [p6| , 

model. 


PGRV corresponds to the calculation of SF F 2 in according to the ref. [29 


Subroutines DF2H8 and DF2D8 calculate F 2 for proton and deuteron within the ref. 
model. 

RANUGL gives the relation between F 2 for deuteron and F 2 for carbon and oxygen; 
F2SFUN is subroutine managing the calculation of F 2 . 

Subroutine R1990 calculates R from ref. [0. 

Fi is calculated in subroutine FISFUN. 

For the calculation of gi in diherent models, the following subroutines can be used: 
PARPOL calculates gi in ref. model; 

SGHAF corresponds the calculation of gi using ref. 
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24] model; 
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GISFUN is subroutine managing the calculation of gi. 

g 2 equals to zero or is calculated in subroutine G2SFUN within the model of ref. 
Quadruple structure functions are calculated in the subroutine B14SF. 

Subroutines FFPRO, FFDEU, FFHE3 and FFGO correspond the calculation of proton, 
deuteron, ^He, carbon and oxygen formfactors. 
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Subroutine FFQUAS calculates formfactors taking into account quasielastic suppression. 
SUPST calculates suppression factors as in ref.[Q; 


PORTN joins together models for SF calculations in different kinematical regions. As a 
result the continuous £t for all SF in whole kinematical region is obtained. 

Patch FITS2 contains subroutines for fit modeling based on the different methods: 

AMNK, GRAM, GAUSS, FI, BASS — on the method of the minimal squared deviations; 

ADIDI, DIVDIF — on the Newton interpolation method; 

GOEFSP — on the cubic spline method; 

MINSTA — on the method of minimal squared deviation by MINUIT (requires GERNLIB). 

Subroutine REMNK2 reads a data hie, calls htting subroutines and stores the parameters 
of constructed hts. The constructed £t can be called by subroutine FITFUN in any place of 
POLRAD. 

Patch INTEGRAT. This hie contains standard integrators DOIFGE, from FORTRAN library 
NAGLIB which are used for the double integration over the square region. Single integration 
is carried out by QUNG8 by the Newton-Gotes method of 8th order, DQG32 by the Gauss 
method or SIMPS by Simpson method. 

Patch POLRAD_ADD corresponds to the calculation of electroweak RG calculation, order 
calculation and approximate RG calculation. It consists on the set of the subroutines (decks), 
that can be divided in three parts: 1) subroutines managing the calculation, 2) subroutines 
that are integrands and 3) auxiliary subroutines. 

1) APPTAI manages approximate and electroweak (only if partonic distributions are 
dependent) calculations; 

TARGWS and GWS manage electroweak corrections calculation when partonic distributions 
are independent; 

AL2LL manages the calculation of corrections contribution; 

2) Here we present the list of subroutines, which calculate the integrands of the corre¬ 
sponding formulae: PEAKl (|2|), PEAK2 (||), UPRE (0|2]), ELU (g|), ELP (||,||), ELQ 
(§§), RA2IPP (|§), RA2ISS (||), RA2ISP (]||), ELUAL2 (0), ELPAL2 (§§, RA2LSS d^), 
RA2LPP (H), RA2FSS (||), RA2FPP (|5|), FXI (1^). 

3) SIGMAB and BOURSG calculates the quantities (H,0) and (|9| JT^) . VERGON dehne 
electroweak coupling constants; SIGALL, VERTS and TT5 calculates one loop ehects: polar¬ 
ization operators (|66D , vertex functions and boxes (eqs.(lO) and (12) of ref. [0) respectively. 
FHOLL and DLAMB are auxiliary function defined in (B.1,B.4,B.6) of ref.[^ VGONEW calcu¬ 
lates born cross section with and without taking into account loop ehects. FFVAPM calculates 
functions (|A.10|) . DISTR calculates partonic distributions. STRFP2 is auxiliary subroutine 
corresponding to the calculation of structure functions. 


3.1.2 Input, Output files 

User can set input parameters for POLRAD 2.0 run in hies POLRAD20.GRA, INPUT.DAT 
and ITDAT*.DAT. 

Input file POLRAD20.CRA contains switches to get necessary case of a calculation. Below 
switches are given with short comments. 

This group of switches corresponds the names of patches, which have to be set on if the 
correspondent part of the correction is calculated: 

polrad - switch on inclusive DIS run of POLRAD ; 
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polrad-add - gives the opportunity to calculate the electroweak RC, order 
effects and QED RC by approximate methods; 
sirad - switch on SIDIS run of POLRAD; 
strffun - calculates SFs (always on); 
fits2 - launches the htting procedure; 
integral - adds subroutines for integrations; 
exact - adds subroutines for exact calculation of RC to DIS. 

The next switches have to be used to choose the way of calculations, type of leptons and 
target polarization: 

approx - gives opportunity to use approximate methods together or instead of 
exact calculation of RC to DIS; 

alpha2ll - calculates the order correction in leading log approximation; 

elect - switches on the electron as an input lepton; 

muons - switches on the muon as an input lepton; 

long - switches on the longitudinal polarization of a nuclear target; 

Iran - switches on the transverse polarization of a nuclear target. 


The next switches are necessary for selection of kinematics within iteration procedure. In 
this case the input data are taken from ITDAT*.DAT. 


iter_pr - launches the iteration procedure for gi or bp, 

minuit - sets the htting of experimental data by MINUIT (only in the case of 
iteration procedure run for gi)-, 

err_prop - calculates error propagation factor in according to eg. ( P-02| ). 
iter_pr_g2 - launches the iteration procedure for g 2 . 


The next switches are necessary for selection of kinematics beyond iteration procedure. 
The only one has to be set on. It should be noted that user can also specify kinematics in the 
beginning of main program body. 


kimnet - switches on the kinematical net over x, y. 
kinsmc - switches on the SMC kinematical set; 
kimhermes - switches on the HERMES kinematical set; 
kin-el42- switches on the E142 kinematical set; 

kimown - switches on the user dehned kinematical set (see hie INPUT.DAT 
description). 

The following switches allow to select the type or nuclear target. The only one has to be 
set on. 


targ_h - switches on proton target; 
targ_d - switches on deuteron target; 
targJied - switches on Helium-3 target; 
targ_c - switches on carbon target; 
targ_o - switches on oxygen target. 
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The next keys allow to select models for SF, partonic distribution and fragmentation func¬ 
tion calculation (see |A.3| also). They correspond to 


f2nmc-d8 - the model for F 2 from ref. [0^ 
f2comfst - the model for F 2 from ref. |^, ^ 

f2glsch - the model for F 2 and gi from ref. 
f2glgrsv96 - the model for F 2 and gi from ref. 


24 


flqpm - the Callan-Gross relation for Fi = F 2 / 2 X] 
r_eq-0 - R = 0] 

glasym- the calculation of gi when iteration procedure is off and switch f2nmc_d8 
or f2comfst is set on; 

qdstr_gu - the model of ref. |Q for partonic distributions; 
g2_eq_0 - 5'2 = 0 ; 

g2-WW - the model for g 2 from ref. ; 

jfrg_aub - the model for fragmentation function from ref. |p5|| ; 


ffrg_cmb - the model for fragmentation function from ref. [pf 
jjrg_arn - the model for fragmentation function from ref. p7 1. 


The next switches allow to specify the quantities to be calculated: 


born - the only Born DIS cross section is calculated; 

poLasym - polarized and unpolarized parts of the total cross section are calcu¬ 
lated separately; 

qua_asym - quadruple and unpolarized parts of the total cross section are calcu¬ 
lated separately; 

cr_sec - total cross section is calculated for all polarized state; 
onlyin - switch to exclude the contribution of elastic and quasielastic tails to 
RC; 

outpuFa - sets asymmetry as a SIRAD output; 

outpuFr - sets quantity r(z) as a SIRAD output; 

intdy - sets the integration of cross section over y; 

intdz - sets the integration of cross section over z] 

cuts - applies the kinematical cuts for RC in semi-inclusive case. 


The next keys guides the calculation of electroweak effects: 

electroweak - calculates the total radiative correction including electroweak RC 
using partonic distributions; 

ew-onlyqed - calculates the electromagnetic RC using partonic distributions by 

eq-dO); 

ew-onlylep - calculates the RC to leptonic current using partonic distributions 
by eq.(IO); 

eweak- calculates the BORN electroweak RC to SIDIS; 

The last group of keys serves for selection of hnal hadron type in the case of SIDIS. They 
switch on the following particles as registered hadron for SIDIS: 
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proton - proton; 
a^proton - antiproton; 
k_minus - K~] 
k_zeroJ)ar - 
k_plus - K^] 
k_zero - K^] 
pi-minus - 7r“; 
pLplus - 7r+; 
pi-zero - 7r°; 

pi-diff - the difference between 7r+ and vr” production to be the measured ob¬ 
servable for SIDIS. 


Input file INPUT.DAT contains eight lines for DIS and nine lines for SIDIS cases. They 
correspond to lepton momentum, target momentum, lepton and nucleons polarization degrees, 
quadrupolarization degree (for deuteron target). The rest lines are used only for kimown switch: 
number of (x, y) pairs for DIS (triads (x, y, z) for SIDIS) and arrays of these pairs (or triads). 


Input file BRASSM.DAT gives coefficients of Brasse ref.|^ for construction of SF fits in 
resonance region. 


Input file BB1FIT.DAT gives coefficients for construction of quadruple SF hl{x) fit |38[ 


Input file STDLOAl.GRI is used for the calculation by model of polarized partonic distri¬ 
butions of ref. . 

Input-output files ITDAT1.DAT, ITDAT2.DAT, ITDAT3.DAT, ITQUAD.DAT, 
ITASM2.DAT contain input-output information for iteration procedure for the cases of ex¬ 
traction of Qi, gf, bi and g 2 respectively. All of them are organized in the same way. Each 
line if it is not a comment (a symbol (not ’O’) in the first position) gives the information on the 
one kinematical point to be processed: x, y (or —Q"^ if negative), measured asymmetry, last 
step extracted asymmetry, previous step extracted asymmetry, measured error, output factor 
for error recalculation given by (|102| ). 

It should be noted that a kinematical point can be removed from analysis if to comment this 
line typing a symbol in the first position. If the symbol is zero correspondent point is processed 
but does not participate in fitting procedure. 

Output files ALL.DAT, ASM.DAT, TAILS.DAT have a title lines with information about 
version, switches used for the calculation in files POLRAD20.CRA and INPUT.DAT. 

The file ASM.DAT gives for each calculated kinematical point quantities x, lU^, Born, 
observed asymmetries and difference between them. 

The file ALL.DAT gives a technical information for each kinematical point about quality of 
numerical integration of all tails and their parts separately for polarization and unpolarization 
parts of cross section. 

The file TAILS.DAT gives eight terms contributed to radiative correction to asymmetry 


Output files ALLP.DAT, ALLU.DAT are usually used to plot a different output quan¬ 
tities. Each line of these files gives the quantities x, ?/, s, 

^eiu,p^ ^eiu,p^ u,p^ respectively. 

It should be noted when the electroweak effects are not calculated the approximate results 
(if approx is set on) for cr“’^ are printed instead of ct“’ 


UyP 

ew' 
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Output files DATA.DAT, FIT.DAT are intended for the control of fitting quality. First 
file contains base points for the fit and the second contains points of the fitting curve. 

3.2 Some examples 

In this subsection we give three basic examples to illustrate the POLRAD 2.0 run: iteration 
procedure of data processing with MINUIT in the case of helium-3 target; the calculation of 
the radiative correction factor for collider DIS; semi-inclusive RC with and without applying 
kinematical cuts. 


3.2.1 Iteration procedure with random input 

Here we demonstrate the RC procedure run within experiments on DIS with polarized Hehum-3 
target with HERMES beam energy E = 27.5GeV. Values of measured asymmetry as well as 
averaged x and were obtained randomly. The following switches were active: polrad, strffun, 
integrat, fits2, exact, elect, long, iter_pr, minuit, targJie3, f2nmc-d8, g2-eq-0, poLasym. 

The output hie ITDAT3 after one step of iteration procedure is given in Appendix p.l| . 
Fig.^ shows these results together with constructed ht of neutron asymmetry (see Appendix 

B- 


3.2.2 Radiative correction in experiments at collider 

Here we give results for unpolarized and polarized radiative correction factors 


du = 


< + af + a% 


(Jn 


dp = 


ag + a? + af 

<^s 


( 112 ) 


within kinematics close to future polarized experiments at HERA collider (lepton and proton 
beam energies equal to 27.5 GeV and 810 GeV respectively). The following switches were 
active: polrad, strffun, integrat, polrad^add, fits2, alpha2ll, exact, elect, long, kin^net, targJi, 
f2glgrsv96, g2-eq-0, poLasym, onlyin, electroweak. The quantities (|112|) are shown in hg.|. The 
output hie ASM.DAT is given in Appendix p.2| (the only x = .001 case is kept). 


3.2.3 Semi-inclnsive radiative correction with and withont cnts 


Here we illustrate the run of the code for semi-inclusive DIS. The following switches were active: 
sirad, strffun, integrat, elect, kin_own, targ_h, qdstr_gu, splineff, jfrg_aub, outfuma, intdz, cuts, 
pi-diff. 

The hg.|^ shows Born and observed asymmetry as well as relative correction 

d = {A°^^ — with and without taking into account kinematical cuts. The output 

hies ASM.DAT with and without applying of experimental cuts are given in Appendices p.3| 
and p.4 


4 Tests and implementation of POLRAD 

POLRAD passed a number of both analytical and numerical tests. It was shown that the 
combinations of coefficients 0(r) coincide with the corresponding combinations of coefficients 
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from ref. |^. Also with the help of the algebraic programming system REDUCE 3.5 we proved 
that the formulae for electroweak inelastic unpolarized radiative correction are the same as 
in Appendix C of ref.[^]. The spin-independent part of POLRAD was tested numerically in 


comparison with FORTRAN codes TERAD86, FERRAD35 and HECTOR and revealed an 
excellent agreement practically in all kinematical regions when all input models and parameters 
were the same. The spin-dependent part of POLRAD was compared with the program of 
Kukhto and Shumeiko |]I| and with E143 radiative correction program kindly placed in our 
disposal by Linda Stuart . We found the satisfactory agreement between all three programs 
when choosing the same models and parameters again. It was impossible to test the part of 
POLRAD corresponding to the calculation of spin one particle quadruple polarization due to 
the total absence of the results in this region. 

Also POLRAD is self-tested program: the part corresponding to approximate calculation 
(see Sect. |2.1.3| ) have excellent agreement with the exact results; the QED part of electroweak 
correction (Sect. |2.2|) coincides with the calculations by the corresponding exact formulae for 
the QED lowest order correction (P!2|); under the simple modihcations of the semi-inclusive 
formulae, the corresponding results coincide with the inclusive ones. 

Now POLRAD is used as the basic and official program for the procedure of radiative 
corrections in SMC (CERN) and HERMES (DESY) and together with the above mentioned 
E143 program in SLAG experiments with polarized particles. 


Appendix A Structure functions 


A.l Inclusive structure functions 


For hadronic tensor we use 






' M ■ *- Jp - 


- (g^„ + 


where 


(A.l) 




3(gh)^ - 

M2 




{p^,Vu + Vi^Pu ) qy 

M2 


, qt^qu ~ f ^ , pqivQ) 

9t.u = 9t.u +Pt^W = + Q2 


(A.2) 


The quantities Aj are dehned as some combinations of SF and formfactors, and the exact 
expressions in dependence of the variant of calculation are considered below in the correspondent 
sections. The dependence of the hadronic tensor on pq and polarization degrees is included in 
Aj too. 

We note, that the hadronic tensor for spin-1/2 particle is derived from ( |A.lj) by putting 
Qat = 0. The formulae for the covariant representation of polarization vector are also valid in 
this case. The hadronic tensor for scalar particles is derived by P/v = Qn = 0. 

Dehning e = /pq we have for various A, in the case of IRT 


29 


























— F -L QiV h 

^1 — ri + Oi, 


% = eiF2 + Sfb2), 


^3 = Pn e{gi + g2), = Pn 5'2, 

^5 = ^6^ bi, % = ^e^b2/3 + b, + bA), 

^7 = ( ^2/3 — &3 ) ; ^8 = ( ^2/3 — &4 ) • 


(A.3) 


Definitions of SF Fj,and bt are the same as in ref. P^ . 

Explicit form of the expression for elastic nuclear formfactors depends on a target spin. For 
deuteron we have 


= IVAPimi + Va) + VaQn), 

%^ = iF^ + hAFl + lv\F^) 

+ ^ {vaFI + F, + F, - F^)F,) , 

55^' = -^(l + 77A)F^(fF, + F,), 

= ^FX2F^ -Fc- ^F,), (A.4) 

C^eZ _ Qn p2 

24 m’ 

^6^ = §f + T4i7(f ) , 

= —^^VAFmiFjn + 2Fq). 


Here = t/AM\ = (Q^ + ReiT)/AM\, and Fc, Fm, Fg - charge, quadruple formfactors of 
deuteron. 

We also give the expressions for the case of arbitrary spin-1/2 nuclei 


C>eZ _ 

Oi — 


Z^VaG 


2 

m’ 


C>e/ _ 
Oo - 




Gl + VaGI 

1 + hA 




C^e/ _ 


PnZ^ 


\Sn = 


GmGe, 


Pi,z^^ G,-Gr, 

^4 — . 1 I _ 1 


1 + hA 


(A.5) 


and for scalar nucleus 

= Z^F^, (A.6) 

where Z is the nucleus charge. All but indicated SF must be set equal to zero. 

The quantities can be obtained in the terms of quasielastic response functions, which 
have a form of peak for uj = Q^/2M. The fact is normally used for construction of the peak 
type approximation. All quantities at response functions are estimated in peak, and subsequent 
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integration of response functions leads to results in terms of suppression factors Se , m,em (or of 
sum rules for electron-nucleus scattering [^). Here we give explicit formulae for ^He target: 

^1 = Vi^^n + 2 /ip)S'M, 


C>9 _ 
^2 — 


+ {el + 2el)S 


l + t] 


^3 2 {Pnen^^n “ 1 “ “^PpCpfip^ SeM i 


C>9 _ 

^4 — 


Pn (-PnCnP'n “1“ ‘2PpCpjJjp') SEM {Pnf^n P ‘2PpH^ SM 

4 1 + y] 


(A.7) 


where rj is rjA for nucleon. Pp and Pn are effective proton and neutron polarization in ^He and 
ep,n, l^p,n are electric and magnetic formfactors of proton and neutron. 


A.2 Electroweak and semi-inclusive structure functions 

Let us define the fermion vertexes and boson propagators. We introduce 

vj = —Cf, aj = 0 , 


_ 7 / ‘^Swef 
7 2s-[ArC\Y 



4 


2sy\rCy/ 


W _ 

f - 


1 

2^/2sw 


(A. 8 ) 


where e/ and — electric charge and the third component of the fermions weak isospin, sw 
and cw — Weinberg angle sine and cosine respectively. Hence, the fermion vertexes takes the 
form 

-ie-ip{v)- (A.9) 
The next couple constants combinations are contained in observables 


\y^ = 2{v'fvj + a'fo’f), 
TDij _ \eij p \eij 


= 2 c/(n}a} + a)v)) , 

p*i _ \eij p \eij 
ri_ — , 


F'Ax, Q^) = x"* ^ [xfxfM(x, Q^) + PkA™i/<-I(i.(J2) 
<? 

F'J{x, Q^) = x"^ ^ [Xf xff*\x, 0 =) + PxK-xff Hx,Q'‘) 


(A.IO) 


where Cf = 1 (— 1 ) for fermions (antifermions), y = Q‘^/{Q‘^ + Ml) and = ( 0 , 1 , 2 ) for ij = 77 , 
■yZ or Zy, ZZ. In QCD-improved model the parton distributions fg'^\x,Q‘^) depend on Q"^: 

ft^{x,Q^) = fP{x,Q^) ± fg^{x,Q^), (A. 11) 
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where and /h(x, — densities of type q partons with helicities, (anti)parallel to 

nucleon helicity respectively. 

The quantities 

S+(-)(x,z) =J2el[f+{x)±f-{x)]Df{z) (A.12) 

Q 

is the ordinary for QPM combination of the distribution functions f^^~\x) for the quark of 
flower q polarized (anti)parallel to the nucleon polarization, and of the fragmentation functions 
Dg{z) of the quark q into the hadron H, Cq being the quark charge in units of elementary 
charge. 


A.3 Fits and models for structure functions 


RC calculation requires hts or models for SF, elastic formfactors, quasielastic suppression fac¬ 
tors, fragmentation functions and partonic densities to be known in a whole region of varying 
of integration variables. 

POLRAD 2.0 gives the opportunity to choose between three models of spin-average and 
spin-dependent proton and neutron (deuteron) SF and gi’'^{x,Q‘^). 

For the explanation of the hrst model for unpolarized SF see hg.^. In the small region 
parametrizations of the ref. (in resonance region) and ref. [^] are used, and for all the rest 
kinematics 15-parameter NA-47 £t is adopted. The advantage of this model is the imple¬ 
mentation of the modern experimental data in small x regions. For the calculation of R{x, Q‘^) 
and F 2 /F 2 the hts shown on the corresponding grafs are used. 

The second model is based on the fit of ||^ with the modern parameters and has the simple 
analytical form, same for the all kinematical region. However, it does not give good description 
of the modern data when x > 0.01. In this model for R we adopt the Whitlow £t 


m 


For both models gi{x,Q‘^) = Fi{x,Q‘^)Ai{x,Q‘^) and Ai(x, could be obtained either 
from the iteration procedure data or from the asymmetry fit (|104|J106D . 

In the third model partonic distrbutions (with [^, ^ or without ^ taking into account 
dependence) are used for the construction of F 2 ’'^{x,Q‘^) and gi{x,Q‘^). This £t is usually 
used in experiments at collider. 

F.^{x,Q‘^) for other nuclei is calculated in accordance with ref. 
to be A-independent 


m 
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and R is considered 

for example). Convolution expressions for gf, 


are 


[see review 
obtained via gi{x) and gi{x) P5| . 

For g 2 one can choose two possible variants: simple partonic approximation g 2 = 0 and the 
Wandzura-Wilczek relation 


92{x, Q^) = -gi{x, Q^) + [ —gi{z, Q^). 

J z 


(A.13) 


Quadruple SF 61,62 should be taken into account for deuterons as spin-1 particles. They 
are related by Callan-Gross equality 


62 ( 3 ^) = 2xhi{x) 


(A.14) 


Model of ref. 38 is used for them. 


and conform to a sum rule obtained in ref. 

Deuteron formfactors are calculated in accordance with the model of ref. , which provides 
the right asymptotic behavior. For ^He form-factors we use £t from ref. . Nucleon formfactors 
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are taken from ref.|^^. Charge formfactor for scalar nuclei can be found in ref.[^] (see sect. 
3.8). 

The suppression factors for QRT for DIS on deuteron target are calculated as the same in 
ref.pQ. For other nuclei we used the Fermi gas model 


There are three possible models for the fragmentation functions 1^^ ^ in SIDIS calcu¬ 

lation. 


A.4 Parameterization of neutron spin asymmetry 


The following function taken from Schaefer’s parameterization is used for htting of the 
neutron spin asymmetry: 


where 


and 


A^x) = 


oo -l- 3ai v 9 


ao — 

ai = 




B(o!u,/^ u-i- 25(0;^,/ 3 d -f 1 ) 
3 

2 B(ad, Pd + 1) 

B = 

J U 


f = _t_ 

“ 1 + auoX°‘'^{l — xY' “ 1 + auoaioX°‘-^ (1 — xY ‘ 

fO ^ _ \ _ ^ 1 

J d -1 I _ _ J d 


(A.15) 


(A.16) 


(A.17) 


^ 1 + ttYttuox'^’^Y — xY^ 1 + aYau{)aiQX°^d.{l — xY 

Only parameters 07 , a^o and Oio are htted, and au, ad, Pu and Pd are considered to be constant 

au = 0.588, ad = 1.03, pu = 2.69, Pd = 6.89. (A.18) 


Appendix B Quantities 


ij' 

In this Appendix we give the explicit form for the functions %(r), which are contained in 
the hnal formulae for the radiative tails, i runs from 1 to 8 . This fact corresponds to the 
contributions of eight SF or formfactor combinations, and j runs from 1 to ki which are dehned 


in (|B.2|) . The function %(r) can be found as a sum over k 


^ij {x ) P ) (^ikTijk{x ) , 
k 

which is calculated from max{l,j + li — kp to min{j, Ip, where 

ki = (3, 3,4, 5, 5, 5, 3,4), k = (1,1,1, 2, 3, 3,1, 2), 

and 


^ik — ^ 


{riq, - 1 } 


for fc = 1, i = 1, 2, 3, 7; 

for k = {1, 2}, i = 4, 8 ; 


- 3 (r/g)^ 67 /g,-3} for A: = {1, 2, 3}, i = 5,6. 


(B.l) 


(B.2) 


(B.3) 
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The quantities Tijkir) for /c = 1 take the form 
Tni(r) = -2m^)FiR, 

Ti 2 i{t) = AtFir, 

Ti3i(r) = - 4 F- 2 t^F,, 

T2n(r) = 2{SX-M^Q'^)Fir/M\ 

T22i(r) = ( 2m^SpF2. + SpS^F,^ + 2{S^ - 2 MH)Fjr - rS^^F ^) /2M^, 

T23i(r) = ( AM^F + (4m2 + 2M^t^ - S,T)Fd - S'pFi+ ) /2M\ 

Tzii{r) = im k2i - ^v)Fir, 

n-n (r) = (,,K(4m^Fi - + 2 F/, - Q'^fL + QlF^ ) + 

+A;2^(-8m^F^'' + Am'^F^^ + 2 rjq tF^) - A r]q + A tFjr), 

T33i(t) = t(F|+- J2i - 2fi) + 

2 tFJ + im^Fl" + %Fjl + - CFg), 


UuM = ^§^(2Ff+ Fg-Fg), 

AttiPt 

UiiM = -p^(S.«fe-2?pQUF,H, 

F42i(t) = ^(4m2(2 Fj - F2+- SpFj + SpF|+) - 2 rS.F,, 

-8 & tF,„ + 2 ( 2S, - S, ) F«, + S, ( Q'‘Fl - Q^F-l)), 
Tm{T) = ^^(2Fi(2m2-TSp) + 2fefTFj + 6F«« 

+ ( 0" - rSp ) F|_ -{ Ql-rS, ) F|+), 

T44i(t) = ^^(2Fi-F|_+F|+), 

Trii (r) = -2 ( + 477^2 ^2 rjki T]k2 ) Fjr, 

T721 (r) = -2(3 T]F{ 2 m^F^_ - rjF tF^ + m Fi+) + 6 rjq F^j^ + tFjr ), 
T73i(r) = 2F + T^Fd + 6{7]FF^^ + r]qTF;i-Am‘^Fr), 


(B.4) 
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^811 (r) = {S rik2 + X r]ki ) Fjr, 

Ts2i (t) = {r]K, ( m^F2- - rSpFd ) + m Fjr + m^SpF^_ 

+ {Sr]ki+Xr]k2)F,+ + S,F^R), 

Ts3i{t) = ^{8m^FS-vlCF,^-vqTFd-S.TF2-SpF^^). 

For i = 5 and i = 6 we have 

T{5,e}jiiT) = T{i^2}ji{T)- (B.5) 

The quantities Tijkir) for fc = 2, 3 are calculated as 

T«t(T) = {f^i ^ ff, ^ fS!,, F’„ ^ FSl) + (B,6) 

The second term appears only for z = 5, 6 and k = 2: 

T 

Qik T 


The substitution in the hrst term of (|B. 6 D has to be applied for all F contained in Tjj_ifc_i(r). 
The quantities F with an upper index are obtained in terms of F without the index: 


2F. 


{iM _ 


2 + 




2F 




AFr 




2 + 


AF. 




2 - 


— ( 2 -Fi++ rF 2 _)s{^,^} + F 2 +r{^_^}, 

= {2,Fd + F2+)ts[^^p} + F 2 _r{g^^}, 

(2-^1+ T ^-^2—) (^»7'S{^,r;} T T -^2+T T 

+4(2F" + FdT‘^)SpS^^^pj, 

= {2,Fd + -F2+)(^r;S{5,77} + + -^2-+ T 


+4rFi+s^S{5,^}, (B.8) 

2Fi\ = {AF + T^Fd)sp + F^+rp, 

AF^l = 2{AF + T^Fd)rpSp + F,+ {rl + T^sl)+A{2F,-TF)sl), 

2F^ = F{rp - TSp) + 2FiSp, 

AF^^ = F{rp-rSpf + AF,{rp-TSp)+AFusl. 

The quantities 

~ = '^(^{5,??} ~ + 2c{5,?7} (B-9) 
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are combinations of coefficients of polarization vectors ^ and rj expansion over basis (see 
sect - P .1.11) 

We note that the scalar products from (|B.5|) and (^) are also calculated in terms of the coeffi¬ 
cients 

VQ = -Q^(ar? - br^) + 5'xC^, r]lC = + 4m^)(a^ -h 6,,) SpC^, 

2r]ki = r]K. + rjq, 2 rik 2 = r]K. — rjq, 
k2^ — Qrn^i ffi 2?71 Wcg, = ScL^ Xb^ 2M c^, 

= 2'm?{a^ap 6 ^ 6 ^) 2M‘^c^Cp + Q‘^{a^bp + b^a^) + S{a^Cp + c^a^) + X{h^Cp - 1 - c^b^). 

(B.ll) 

The following equalities define the functions F\ 


F = 


Q ’ 


Fir — F 2 + — QlaFd, 

F, = r-‘(Cj-‘'"(T) - cP'\t)) Fi+ = cP'^(t)+cP'\t), 

F 2 + = B^{t)CP''\t) - Bi(t)CP''\t), Fj- = S2(t)C2-=‘''"(t) + B,{t)cP''‘{t) 


(B.12) 


F, = 


= (3BiW-AeCi(T)), 


where 


5i,2(r) = -- ( Aqt ± Sp{S^T + 2Q2)), 
C'i(r) = {Sr + Q2)2 + 4 ^ 2 (g2 + 


(B.13) 


C2{r) = {Xt - g 2)2 tS^ - 

We note that F^ has a uncertainty like 0/0 for r = 0 (inside of integration region). It leads to 
difficulties for numerical integration, so the another form is used also 


Fb = 


Sp(tS, + 2 Q= 


C;-'"(r)Q'^(r)(C;"(T) + C 2 ‘"(T)) 


U/2 


1 / 2 / 


W 2 / 


(B.14) 


Appendix C Quantities T^i’ 


u,p 


The functions contributed to (|56|) are listed here: 




-2Y,Xi + 2yv((| - l)f(r,) - ; + 3 ^) - ® > 


+(f+^//(.ri-3)/h 
3 3 2yip X 


— (1 — x){ -h Y+)Li + —{2yi + Y+)Lyi 

-’2 y^ 2p 
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+ ^(~ + ~ ~ ~ ^ ^+)Ly + 


7^“ = -2F+L + |-(2a;/(r/)-3) + (|- + F+)Li + (2|/i + F+)L^i-4y+L, 

~y{Dl + -Dyl), 

7^? = 2Y_xXL-2Y(l-- + 2fiT])) + 2iY(l--)-—)Li + 

X X xyi 

2(r_(i - i) + ^)L„ + ^(y_ + ^)Li - 5 y_y+Li - |^(2y+ - Yl,i 

X X 2rj Hi Tj 27] yi 

+ {2y + 2 --- )Di + (2- 2.y — - - )Dyi, 

X 2yi7] X 2yi7] 

nl = 2Y_L + Y_{L,-l + Ly,)+y{D,-Dy,), 


(C.l) 


where 


Ly = \n{xf{y)), Li = \n{yi + xyf{y)), Lyi = ln( 


1 - xyfjy). 
yi 


Di = 
L = In 




yi+xyf{yy l-xyf{yy 

(1 — xy 


(1 - xf{7])y{l - xy){yi + xy) 




yi 


and 2/(77) = 1 + \J^h + 1- 


(C.2) 


Appendix D Test Run Output 
D.l Example 1. Output file ITDAT3.DAT 

*************** Heliuiii~3 data ********************** 


010 

.830 

-.03500 

-.04566 

-.03500 

.02000 

020 

.790 

-.05500 

-.06745 

-.05500 

.02100 

035 

.740 

-.07500 

-.08794 

-.07500 

.02300 

050 

.660 

-.10000 

-.11211 

-.10000 

.02700 

070 

.620 

-.08000 

-.09195 

-.08000 

.02800 

095 

.570 

-.12800 

-.13951 

-.12800 

.03000 

115 

.530 

-.14000 

-.15128 

-.14000 

.03200 

140 

.480 

-.16000 

-.17111 

-.16000 

.03400 

180 

.450 

-.12500 

-.13620 

-.12500 

.03500 

230 

.425 

-.17000 

-.18141 

-.17000 

.03600 

285 

.400 

-.14500 

-.15666 

-.14500 

.03800 
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.350 

.360 

-.12500 

-.13701 - 

.12500 

.04000 

.430 

.330 

-.08000 

-.09231 - 

.08000 

.04200 

.530 

.300 

.00800 

-.00438 

.00800 

.04400 

.640 

.270 

.09000 

.07667 

.09000 

.04600 

.745 

.240 

.22000 

.20203 

.22000 

.04800 

.850 

.220 

.48800 

.45913 

.48800 

.05000 

Next two 

lines 

are additional data 

for fit 

constructing 

* 1.000 

0.000 

1.00000 

1.00000 1 

.00000 

.00000 

* .000 

.000 

.00000 

.00000 

.00000 

.00000 


Am - measured spin asymmetry, 

A1 - last step spin asymmetry, 

Ap - previous step spin asymmetry. 

Err- experimental error. 

D.2 Example 2. Output file ASM.DAT 

program polrad20 version from 10.04.1997 

the file gives born asymmetry, observed asymmetry 
and radiative correstion 

the following switches are active 

polrad strffun integrat polrad_add fits2 alpha211 exact elect long 
kin_net targ_h f2glgrsv96 glasym g2_eq_0 pol_asym onlyin 
electroweak 

leptons are electrons 
target is proton 

target is longitudinally polarized 
bmom = 27.5 

tmom = 830.0 


pi = 

1.00 pn 

= 1.00 

qn = 

0.00 


is in 

7o 





X 

w2 

q2 

a(born) 

a(obs) 

del(%) 

0.001 

913.0 

0.9 

0.378 

0.334 

-0.043 

0.001 

1369.0 

1.4 

0.378 

0.337 

-0.040 

0.001 

1825.1 

1.8 

0.378 

0.339 

-0.038 

0.001 

2737.1 

2.7 

0.378 

0.342 

-0.035 

0.001 

3649.2 

3.7 

0.378 

0.345 

-0.032 

0.001 

4561.3 

4.6 

0.378 

0.347 

-0.030 

0.001 

5473.4 

5.5 

0.378 

0.349 

-0.029 

0.001 

6385.5 

6.4 

0.378 

0.349 

-0.028 

0.001 

7297.6 

7.3 

0.378 

0.350 

-0.027 

0.001 

8209.7 

8.2 

0.378 

0.351 

-0.027 
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0.001 9121.8 
0.001 11402.0 
0.001 13682.2 
0.001 15962.4 
0.001 18242.6 
0.001 27363.5 
0.001 36484.4 
0.001 45605.2 
0.001 54726.1 
0.001 59286.6 
0.001 63847.0 
0.001 68407.4 
0.001 72967.9 
0.001 75248.1 
0.001 77528.3 
0.001 79808.5 
0.001 82088.7 
0.001 83000.8 
0.001 83912.9 
0.001 84825.0 
0.001 85737.1 
0.001 86649.2 
0.001 87561.3 
0.001 88473.3 
0.001 89385.4 
0.001 90297.5 


9 

.1 

0 

.378 

11 

.4 

0 

.378 

13 

.7 

0 

.378 

16 

.0 

0 

.378 

18 

.3 

0 

.378 

27 

.4 

0 

.378 

36 

.5 

0 

.378 

45 

.7 

0 

.378 

54 

.8 

0 

.378 

59 

.3 

0 

.378 

63 

.9 

0 

.378 

68 

.5 

0 

.378 

73 

.0 

0 

.378 

75 

.3 

0 

.378 

77 

.6 

0 

.378 

79 

.9 

0 

.378 

82 

.2 

0 

.378 

83 

. 1 

0 

.378 

84 

.0 

0 

.378 

84 

.9 

0 

.378 

85 

.8 

0 

.378 

86 

.7 

0 

.378 

00 

.6 

0 

.378 

00 

00 

.6 

0 

.378 

00 

.5 

0 

.378 

90, 

.4 

0 

.378 


0 

351 

-0 

027 

0 

352 

-0 

026 

0 

352 

-0 

025 

0 

353 

-0 

025 

0 

353 

-0 

024 

0 

355 

-0 

023 

0 

355 

-0 

022 

0 

355 

-0 

022 

0 

356 

-0 

021 

0 

355 

-0 

023 

0 

353 

-0 

025 

0 

350 

-0 

027 

0 

346 

-0 

032 

0 

343 

-0 

035 

0 

341 

-0 

037 

0 

336 

-0 

041 

0 

336 

-0 

041 

0 

324 

-0 

054 

0 

323 

-0 

055 

0 

320 

-0 

058 

0 

310 

-0 

067 

0 

298 

-0 

079 

0 

291 

-0 

087 

0 

278 

-0 

099 

0 

261 

-0 

116 

0 

235 

-0 

142 


D.3 Example 3. Output file ASM.DAT (with CUTS) 

program polrad20 version from 10.04.1997 

the file gives born asymmetry, observed asymmetry 
and radiative correstion 

the following switches are active 

sirad strffun integrat elect long kin_net targ_h f2glgrsv96 glasym 
qdstr_gu g2_eq_0 ffrg_aub pol_asym onlyin outfun_a intdz cuts eweak 
pi.diff 

leptons are electrons 
target is proton 

target is longitudinally polarized 
bmom = 27.5 

tmom = . 0 

pi = 1.00 pn = 1.00 qn = .00 
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a is in % 


X 

y 

q**2 

A1 

rc A1 

born 

meas 

115 

.421 

2.490 

58.912 

-2.828 

30.237 

29.382 

183 

.421 

3.980 

61.646 

-1.666 

31.289 

30.768 

252 

.421 

5.470 

63.000 

-.986 

31.780 

31.467 

320 

.421 

6.960 

63.683 

-.574 

32.015 

31.831 

389 

.421 

8.449 

63.999 

-.322 

32.116 

32.013 

457 

.421 

9.939 

64.005 

-.198 

32.103 

32.040 

526 

.421 

11.429 

63.979 

-.171 

32.106 

32.051 

594 

.421 

12.919 

63.938 

-.298 

32.150 

32.054 


D.4 Example 3. Output file ASM.DAT (without CUTS) 

program polrad20 version from 10.04.1997 


the file gives born asymmetry, observed asymmetry 
and radiative correstion 


the following switches are active 

sirad strffun integrat elect long kin_net targ_h f2glgrsv96 glasym 
qdstr_gu g2_eq_0 ffrg_aub pol_asym onlyin outfun_a intdz eweak 
pi.diff 


leptons are electrons 
target is proton 

target is longitudinally polarized 


bmom = 

27.5 






tmom = 

.0 






pi = 

1.00 

pn = 1, 

.00 qn = 

.00 



a is in 

7o 






X 

y 

q**2 

A1 

rc A1 

born 

meas 

.115 

.421 

2.490 

57.759 

-4.730 

30.237 

28.807 

.183 

.421 

3.980 

60.751 

-3.093 

31.289 

30.321 

.252 

.421 

5.470 

62.364 

-1.985 

31.780 

31.149 

.320 

.421 

6.960 

63.256 

-1.241 

32.015 

31.618 

.389 

.421 

8.449 

63.742 

-.722 

32.116 

31.885 

.457 

.421 

9.939 

63.908 

-.349 

32.103 

31.991 

.526 

.421 

11.429 

64.045 

-.069 

32.106 

32.084 

.594 

.421 

12.919 

64.212 

.129 

32.150 

32.191 

.663 

.421 

14.409 

64.189 

.303 

32.106 

32.204 

.731 

.421 

15.899 

64.364 

.481 

32.160 

32.315 
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FIGURE CAPTIONS 


1. The double differential cross section for inclusive lepton-nuclei scattering sketched at a 
certain value of as function of u (arbitrary scale). The three basic channels are the 
elastic (I), the quasielastic (II), and the inelastic (III) one. 

2. The limits of integration for variables a) R and r (|T^; b) ^ and t 

m- 

3. Scheme of program POLRAD. The subroutine APPTAI, AL2LL and TARGWS are pre¬ 
sented on a separate hgures. 


(|3^^^ ; b) ti and t 2 


4. 

5. 

6 . 

7. 

8 . 


Scheme of program SIRAD. 

Scheme of program AL2LL. 

Scheme of program APPTAI. 

Scheme of program TARGWS. 

The results of iteration procedure for spin asymmetry (see sect. 3.2.11) along with the £t 
constructed. 


9. The unpolarization a) and polarization b) radiative correction factors defined by eqn.( 112 ). 
The curves 1,2 and 3 correspond to different values of a; = 0.001, 0.01, 0.1 respectively. 


10. The born (dash) and observed (solid) SIDIS asymmetries (a) and relative correction with 
(dash) and without (solid) kinematical cuts apllied (b). 


11. The spin-independent proton and neutron SF. The citations denote that fit from the 
cited references is used in the range. The fits are extrapolated into hatched region. The 
procedure of joining together of two-dimensional surfaces gives continuous fits in the whole 
kinematical region. 
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